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Abstract 

High accuracy helium wave functions based on exponentials with random coefficients are trans- 
formed into momentum space. The utility of the wave functions is demonstrated through cal- 
culation of the expectation value of various operators needed to evaluate relativistic and QED 
corrections. 
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I. INTRODUCTION 



In recent years a basis set consisting of random exponentials has been used with increasing 
frequency to carry out calculations in helium □ no . For singlet and triplet S-states it 
can be written as 



0(f 1; f 2 ) = c i [ e - Qiri - ftr2 - 7 ' ri2 ± e - ain - piri -* na ], (i) 

i 

where «j, and % are parameters that are chosen randomly in certain ranges, and the 

n 

spin wavefunctions are understood. Through careful choice of those ranges Korobov [l( has 
been able to obtain a ground state energy of 

E = -2.903 724 377 034 119 598 311 159 a.u.. (2) 

While even higher accuracies are possible with basis sets that incorporate known nonanalytic 
behaviors of the wave function 4J, the simple form of the above wave function makes the 
evaluation of higher order corrections coming from relativistic and QED corrections relatively 
straightforward. In this note we will present calculations based on the momentum space form 
of Eq. 1, where the wave function is defined through 

<P(PUP2) = J ^^e-^e-^^(n,r 2 ). (3) 

The Fourier transform is evaluated by first noting that 

^3 e -otr 1 -f3r 2 -7r 12 



e -ari-/3r 2 -7ri2 



d 3 f d 3 qi e** 1 * 1 f d 3 q e i?f?12 f d 3 q 2 e^ 2 



(4) 



dad(3d-f J 2tt 2 q x 2 + a 2 J 2vr 2 q 2 + 7 2 J 2vr 2 q 2 2 + (3 2 ' 

To simplify the following discussion, we Fourier transform only the above expression, and 
define it as 0j(pi,p 2 ), with the generalization to Eq. 1 being clear. We see that 

1 d 3 f 3 1 1 1 

&(pi,£0 = ~^3 dad/3dl J d Qfa _ £)2 + a 2 + 7 2 + ^2 +(3 2- ( 5 ) 



The integral over d 3 q has been carried out analytically in Ref. j^, and leads to 

, ^ 2 d 3 6 

4>i\P\,P2) = o a nr> 7=, (6) 

Ti oaopoj \JX 
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where 



V = l[\P2 + Pi\ 2 + {a + (3) 2 } + /5( 7 2 + a 2 + pi 2 ) + a( 7 2 + (3 2 + p 2 2 ) 
x = [\P2 + Pi I 2 + (a + P) 2 ] \pi 2 + (a + 7) 2 ] [p 2 2 + (/3 + 7) 2 ] - Z/ 2 



9 = arctan(^). (7) 

While this equation for <pi(pi,p 2 ) seems quite compact, the action of the three derivatives 
leads to a considerably more complicated expression. It simplifies if we introduce the auxil- 
iary parameters 

dx dx dx 

1 da 2 dp 1 3 g?7 

Yi = ^L,Y 2 = ^-,Y 3 = ^. 
da d/3' c?7 

(i 2 x d 2 x <i 2 x 

X 2 3 = m ) -^12 = 1 77^X13 



dpdj dad(3 dadj 

T\ = XiX 2 3 + X 2 X 13 + X 3 X 12 , T 2 = X\ + X 2 + X 3 , T 3 = yX^ + Y 2 Xi 3 + Y 3 Xi 2 

T 4 = Y 1 + Y 2 + Y 3 , T 5 = Xiy 2 y 3 + X 2 Y 1 Y 3 + X 3 Y l Y 2 

Tq = Y\X 2 X 3 + Y 2 X\X 3 + Y 3 X\X 2 . (8) 



In terms of these we find 



0i(Pi,P 2 ) = --[g^pa^ - 5^X2X3] + — + — 



4yT 5 2(x - 3 j / 2 )y 1 y 2 y 3 (40x?/ 3 + 33x 2 ?/ + 15 j / 5 )X 1 X 2 X 3 
L> 3 D 3 8x 3 D 3 

yTi(3y 2 + 5z) 7T X 2 1 



(9) 



4x 2 D 2 D 3 D 1 ' 

where D = x + y 2 and y" = 2(a + (3 + 7). We note the following symmetry properties of the 
momentum space wave function. Because the basic functions that form the wave function, 
x and y, are invariant under the simultaneous replacement p\ <-> p 2 and a (3, for singlet 
states 4>(pi,p 2 ) = 4>(p 2 ,Pi) and for triplet states 4>(pi,p 2 ) = —<P(p 2 ,Pi)- Both symmetries 
also hold when the magnitudes of p\ and p 2 are switched. 

There are various uses for the wave function of helium in momentum space, notably 
application to scattering calculations. While most work has been carried out with Fourier 
transforms of Hartree-Fock wave functions, an approach that includes correlation more com- 
pletely is given in was Fourier transformed and 
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applied to the calculation of a number of helium properties, notably the Compton profile. 
While the wave function used here gives a more accurate energy (odd powers of ri2 were not 
included in the basis set of Ref. [6j]), we do not find appreciably different answers for any 
of the quantities calculated there. In particular, we find the same small difference between 
using fully correlated wave functions and Hartree-Fock wave functions for the Compton 
profile. Rather than pursuing this line of research we instead now discuss applications of 
the momentum space wave function to the calculation of higher order relativistic and QED 
corrections to the energies of the l 1 ^, 2 1 S , and 2 3 Si states of neutral helium. 

These corrections are obtained from a set of operators Oi, with associated energy shifts 
Ei = {(j)\Oi\(f)). The first set of operators was derived by Breit 7], and describe corrections of 
order a 2 a.u.. The equation he used in the derivation has certain difficulties connected with 
negative energy states g but later treatments using the Bethe-Salpeter elation Q, Q, 



111 ] treat negative energy states consistently and allow the systematic treatment of higher 



order corrections. 

Most recent calculations, however, use the technique of effective field theory |12| to derive 
the operators. We note in particular the compact rederivation of the Douglas-Kroll ^j] 
results for contributions to the fine structure of helium in order a 4 a.u. of Ref. Q], and the 
derivation of contributions to the energy of the ground state to the same order in Ref. jl^j ]. 
The idea of effective field theory is to compare free-particle scattering amplitudes in QED 
to an effective nonrelativistic theory, with operators added to the Schrodinger Hamiltonian 
that account perturbatively for the difference of the amplitudes. Once the operators have 
been determined from considering free-particle scattering, they are used as perturbations in 
standard bound-state Rayleigh-Schrodinger perturbation theory to calculate energy shifts. 
In this approach it is natural to work in momentum space, and the operators are then 
Fourier transformed to coordinate space for numerical evaluation. Here, however, because 
we have formed wave functions in momentum space, we avoid this step, and work exclusively 
in momentum space. This has the advantage of simplicity, but the disadvantage of being 
less accurate than coordinate space. We now give a brief rederivation of the Breit operators 
using effective field theory, and illustrate their numerical evaluation in momentum space. 

The connection of the scattering amplitude of two electrons with momentum pi, p 2 to 
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scatter into states with momentum p 3} p^ to an energy shift is 

Ei = j^-r^ J d 3 p 1 d 3 p 2 d 3 p 3 d 3 p^(/)%p 3 ,p4)Oi(p 3 ,p4]Pi,P2)(p(pi (10) 

We work in the center of mass frame, with initial nuclear momentum —p\ — p2, final nuclear 
momentum — p 3 — p 4 . If only electron-electron interactions (ee) are considered, momentum 
conservation allows us to write 

Oi ee {p3,P^Pi,p 2 ) = (2ix) 3 5 3 (p 3 +p A -p 1 -p 2 )Mr(p3;PuP2) (11) 

with the associated energy shift 

Ei ee = J d 3 p 1 d 3 p 2 d 3 p 3 (f)*{p 3 ,p 1 +p2 - p 3 )^i ee (P3;Pi,P2)0(pi,^2)- (12) 

If instead we consider diagrams in which one of the electrons, taken to be electron 1, interacts 
with the nucleus (eN), and electron 2 does not participate, we can write 

Oi eN (p 3 ,p 4 ;puP2) = (2rt) 3 6 3 (p 4 - p^M^^Pi) (13) 

with the energy shift 

Ei eN = j^y 3 J rf 3 p 1 ^ 2 ^ 3 0*(P3,P 2 )M/^(pl 5 ,^)0(p 1 ,p 2 ). (14) 

Diagrams with all three particles interacting have no delta functions, and have to be evalu- 
ated with the 12-dimensional integral of Eqn. El In all cases we note that 3 of the integration 
variables can be carried out trivially, and that the adaptive Monte-Carlo program VEGAS 
[lsl can be used to numerically evaluate the integrals, though with far less accuracy than 
available from coordinate space techniques. This loss of accuracy is due entirely to the fact 
that the multidimensional integrals have to be carried out numerically: the wave functions 
themselves are quite accurate. In the calculations presented here we use 600 basis functions, 
and the energy eigenvalues are accurate to more than 14 digits for the ground state and 11 
for the excited S states. 

We begin by evaluating the relativistic mass increase (RMI) operator, which we treat as 
an ee diagram. The contribution to the scattering amplitude of this operator is 

M RMI ee (p 3 ;p u p 2 ) = -(2n) 3 ^£^5(p 3 - pi), (15) 
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which gives the energy shift 

E RM i ee = -g^g / d 3 Pl d 3 p 2 <f>*(pi,P2)(pi +pMpuP2)- (16) 

The result is tabulated in the first row of Table 1. 

We next turn to corrections to Coulomb scattering between the electrons. In this case 
the nonrelativistic scattering operator is 

Al'KCH 

M c ee (p 3 ;Pi,P2) = p (17) 

IPs -Pi\ 2 

This corresponds to the coordinate space potential a/\fi — f 2 \. To calculate relativistic 
corrections, we use Dirac spinors to describe scattering and work in Coulomb gauge. We 
introduce the notation DC to refer to the scattering amplitude with exchange of a Coulomb 
photon, and DT for the scattering amplitude with exchange of a transverse photon. The 
DC scattering amplitude can then be Taylor expanded in powers of p/m, with the leading 
corrections given by 

M DC ee (p 3 ;PuP2) = Mo-Cft; ft, - |P Vf ' 2 " b \~? 12 + 

ioi ■ (p 3 x gi) iff 2 ■ (p 4 x p 2 ) 
4m 2 4m 2 

with the understanding that p^ = p\ + p 2 — p$. The first two correction terms are Darwin 
terms, and sum to na/m 2 S 3 (f2 — fi) in coordinate space. The last two are spin-orbit op- 
erators, which do not contribute to the S-states considered here. We tabulate the Darwin 
terms in the second row of Table 1. 

Considering now Coulomb scattering between an electron and the nucleus, which we take 
to have charge Z although only Z = 2 will be considered here, the nonrelativistic limit is 
given by 

Mc eN (p\p 1 ) = -r^r V2 (19) 

\P3 -pi\ 2 

with associated energy shift 

E c eN = -j^y 3 J d 3 Pl d*p 2 d 3 p 3 <P* (p 3 , p 2 ) ^ Z - m ,p 2 ). (20) 

Relativistic corrections are now obtained by introducing a Dirac spinor only for the electron 
(the nucleus is treated here in the infinite mass limit), and we find for exchange of a Coulomb 
photon 

Moc-fe.S) = M c -(p s ,p 1 )[l - + i3l -fj ?l) ] (21) 
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Again only the Darwin term contribute for S-states, but now corresponds to TrZa/m 2 S 3 (fi), 
which we tabulate in the third row of Table 1. 

The effect of transverse photon exchange between the electrons is simplified if we neglect 
retardation, which enters in order a 3 a.u., and in this approximation we have 

7TCt 7TCt 7TQ 

M DT ee (p 3 ; pi,p 2 ) = —r 5~2 (Pi + Pa) ■ (02 + P4) A3 x -B 2 -d x - qa 2 ■ q) 

2ma _ 2nia _ ._ , . 

m z q z m z q z 

where q = \p% — p\ \ . The first two terms, referred to as orbit-orbit terms, or as are 
usually evaluated by Fourier transforming into coordinate space. With the present approach, 
however, they are quite easily treated, and the result presented in the fourth row of Table 
1. The third term is another delta function, and the last terms again vanish for the S-state 
considered here. The effect of transverse photon exchange for eN scattering vanishes in the 
infinite nuclear mass limit used in this work. 

The operators considered so far in this note have been studied for many decades, and have 
all been evaluated with much higher accuracy than presented here [16]. The utility of the 
present approach lies in the fact that operators that enter in higher order, generally derived 
in momentum space, are both fairly complicated when Fourier transformed to coordinate 
space, and in addition need to be evaluated with less accuracy than the terms treated above. 
We illustrate this point with relativistic corrections that contribute in order ma(Za) 6 to 

n 

S-states. A complete set of operators for triplet states has been derived by Pachucki |12| 
using an effective field theory approach, and we consider here the corrections to one-Coulomb 
photon exchange, Eq. 20 of that paper, 

= q 2 64m 4 ^ ~ 3^ 3 X Pl ' ' ^ 4 X P2 '^ 1 ' ° 2 ' + 

\U - Pi) 2 + (Pi - Pi) 2 } + m 2 {p\ +pI+pI + pI)}- (23) 

This expression corresponds to the next term in the p/m expansion in Eq. While Ref. 



|17| treats triplet states, this particular result is also valid for the singlet case. We note 
that the last term vanishes for triplet states because of the symmetry <f>(jpi,p2) = —<P(P2,Pi) 
mentioned earlier, as also noted in Ref. Q]. The resulting energy shift of the 2 3 S*i state is 

E ee = -0.0062(2)ma 6 . (24) 
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While again not of high accuracy, we note the extremely simple nature of the coding, which 
is almost identical to the program that evaluates the Darwin term. This contrasts with the 
more complicated coordinate space calculation, where numerous derivatives must be applied 
to wave functions, leading to a much lengthier expression. As with the ma 4 corrections, much 
higher accuracy is available from working in coordinate space, with -0.006 344 7 ma 6 the 
known result 18| . However, we note the momentum space accuracy corresponds to 3.7 kHz, 
to be compared with the experimental accuracy of 60 kHz. 

While the formula for V\ given above is valid for singlet states, it gives a linearly divergent 
result in that case. It is quite simple, however, to regulate this divergence in momentum 
space, where one simply imposes the cutoff \pi\ < A. An application of this momentum space 
regulator to the case of ground state positronium hyperfme splitting can be found in [20(. In 
Table II we show results for the expectation value of V\ for the ground state of helium with 
different cutoffs A, with the linear dependence on A clearly visible. When combined with 
other linearly divergent terms in a complete calculation a A independent result will obtain 
in the limit of large A. By improving the accuracy found in this part of the calculation 
this procedure can be used to check the results of Ref. [^] without explicitly canceling the 
divergences: work on this problem is in progress. 

II. CONCLUSIONS 

We have presented the formula for the momentum space form of a powerful basis set for 
helium. While it has the potential for proving useful for scattering calculations on helium, 
we have concentrated on evaluating expectation values of operators that give relativistic and 
QED corrections to energy levels. Because these operators are derived in momentum space, 
this allows one to work entirely in momentum space. The next step in this research is the 
extension to states with nonvanishing angular momentum. The most important application 
we have in mind is to the fine structure of helium P states, where recent high-precision mea- 



surements by Hessels and collaborators 



2l| have the potential of allowing a determination 



of the fine structure constant a to a precision of 4 ppb. Unfortunately, the present state of 



theory is unclear, where the most complete calcu 
the fine structure interval v$\ measured in Ref. 



ation by Drake [22], while consistent with 



3, is i 



inconsistent with measurements of 



the interval V\2 j^, jsJ). This inconsistency has also been noted in Ref. j^. It is possible 



S 



that the relative simplicity of the method developed here can shed light on this situation. 
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TABLE I: Expectation of operators for n=l and 2 S states of helium in units of a 2 a.u.. The 
notation H2 in E[)T ee indicates the delta function is not included in the result. 



Operator 


l 1 ^ 


2 1 S 


2 3 5i 




-13.5212(3) 


-10.27959(5) 


-10.45887(4) 


E DC ee 


0.3346(3) 


0.02718(8) 


0.0 


E DC eN 


5.6879(2) 


4.1139(2) 


4.1479(2) 


EDT ee {H2) 


-0.1393(2) 


-0.00922(1) 


-0.00157(7) 



TABLE II: Expectation value of V\ for the ground state of helium with the regulator \pi\ < A for 
different values of A. Units a 4 a.u. for (V\) and ma for A. 

A 100 200 300 

(Vl) -26.9(1) -59.4(4) -90.9(6) 
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